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The question of the variation of (1) the direct and (2) the diffused 
radiation is next discussed. On comparing the curve representing the 
chemical intensity of diffused light with the curve of solar obscuration, it 
is found that the rate of diminution in chemical action exerted by the dif¬ 
fused light is up to a certain point greater than corresponds to the portion 
of eclipsed sun, whilst from this point up to totality the rate of diminution 
becomes less than corresponds to the progress of the eclipse. The same 
rapid diminution in the chemical action of the diffused daylight during the 
early periods of the eclipse was also observed at Jamkandi; it is doubtless 
due to the dark body of the moon cutting off the light from the brightly 
illuminated portion of sky lying round the solar disk. 

The results of the observations at Catania are then compared with those 
made at Moita, near Lisbon, and communicated to the Society in 1870. 
This comparison shows a striking coincidence between the two sets of ob¬ 
servations. In each case it is seen that the relation between solar altitude 
and total chemical intensity is represented by a straight line, although the 
Catania observations slightly exceed, by a constant difference, those made 
at Moita in conformity with the slight difference in latitude, and with the 
fact that the former determinations were made at a greater elevation above 
the sea-level. 

The Catania observations further confirm the fact which we formerly 
announced, that for altitudes below 50° the amount of chemical action 
effected in the plane of the horizon by diffused daylight is greater than that 
exerted by direct radiation, and also that at altitudes below 10° direct sun¬ 
light is almost completely robbed of its chemically active rays. 


VII. a On the Calculation of Euler* s Constant.” By J. W. L. 
Glaisher, B.A., F.R.A.S. Communicated by James Glaisher, 
F.R.S. Received June 6, 1871. 

The main object of the present communication is to correct some inac¬ 
curacies both of reasoning and calculation contained in two papers by Mr. 
Shanks, viz. Cf On the Extension of the Value of the Base of Napier’s 
Logarithms ; of the Napierian Logarithms of 2, 3, 5, and JO ; and of the 
Modulus of Briggs on the common System of Logarithms; all to 205 
places of Decimals,” in the Proc. Roy. Soc. vol. vi. p. 397; and "On 
the Calculation of the Numerical Value of Euler’s Constant,” in the Proc. 
Roy. Soc. vol. xv. p. 429 (1867). 

For the calculation of the constant Mr. Shanks has used (as, indeed, 
has every calculator who has computed the value of the constant during 
the present century) the semiconvergent series 


y 


1 i 

2 3 


, 1 1 1 , 

+ s- loga; “^ + 


2x 2 4x A * ‘ * 3 


0 ) 


y being the constant, and B,, B s , B„ . , . Bernoulli’s numbers. 
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Mr. Shanks obtains the value of y from this formula by making x= 10, 
20, 50, 100, 200, 500, and 1000, and remarks, as a curious coincidence, 
that the number of decimal places obtained from x being made equal to 
10, 20, 50, and 100 is nearly proportional to ^/lO, |/20, %/50, and ,3/100. 
On p. 432 he gives a Table of the number of decimal places obtainable 
from the formula when x has the values 2, 5, 10, 20, 50, 100, 200, 500, 
and 1000, and from it draws the inference that “we may fairly infer that 
when n is increased in a geometrical ratio, the corresponding number of 
decimals obtained in the value of E increases only in something like an 
arithmetical one, and that probably from 50,000 to 100,000 terms in the 
Harmonic Progression would require to be summed in order to obtain 100 
places of decimals in the value of E, Euler’s constant.” 

Algebraically, of course, y is independent of x in the formula (i), but 
arithmetically, since the series ultimately becomes divergent, the value of y 
is so far dependent on x that for a given value of x the series will only 
afford a certain number of decimal places. The number of decimal places 
directly obtainable is equal to the number of ciphers which precede the first 
significant figure in the value of the numerically least term of the series. 

The nth term (considering only the terms after ^ in (i), so that the 

nth term is the term involving B„) 

- B n 

2nx 2n ' 

which, since 


B _ 2(l,2.3...2«y . 1 _1 

(2rc) 2n \ n V ,1 3 2 ' 1 

is very nearly equal to 

2{1.2.3 . . . (2n~~ 1)} # 

(27rx) 2rl 

so that the ratio of the nth to the (n—l)th term 



__ (2n—2)(2n— 1)__ n* _ 3 \ 

4ttV k 2 x\ 2 n) 

very nearly. 

Let m be the greatest integer contained in xtt so that ^7r=m~f-/, f 
being a proper fraction, then the ratio of the mth to the (m— l)th term 


nr 


o-&) 


(m+f) 

-iJ/+3 

2m ' 

which is always less than unity. 

The ratio of the (m+l)th to the 'mth term is found in a similar manner 
to be 

i+ s(H- 


2 s 2 
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which is greater or less than unity according as/is less or greater than 
The ratio of the (m-f-2)th to the (m+ 1) th term 



which is always greater than unity. 

Thus the mth or (m-f- 1)th term is the least according as / is less or 
greater than i. 

When m is large the value of the mth term is very nearly 


1 1.2.3_2m 

m 


(2tt x) 2m 

which, by use of the theorem 

1.2.3...®= ^(2rr.v)«fe-*^l+ 2-), 


becomes 


1 . 2 a/ 7r(2m) 2 'V 


sj m Qhzx) 2 


( 1 + i) 

“VI* fe) 2 "‘( 1 + i)- 


(ii) 


The number which forms the negative characteristic of the common 
logarithm of this quantity (the mantissa being made positive), reduced by 
unity, denotes the number of ciphers which precede the first significant 
figure in its value; so that if the mantissa be not made positive, the cha¬ 
racteristic indicates the number of decimals directly obtainable from the 
series. 

The logarithm of (ii) with its sign changed, after some expansions and 
reductions, becomes 

2)XirX + h°S 10* — log 10 2 — /j (J— + 

jL \Zxk xj r Z4X7T/ 


fi being the modulus 43429448 .... Neglecting the last three terms, 
which can only very rarely lead to an error of a unit, we obtain as a result 
that the number of decimal places which the formula (i) will afford directly 
for a given value of x is equal to the greatest integer contained in 

2pir%+~\og 10 %— log 10 2.(iii) 

The corresponding value of the least term is, of course, easily obtained 
from (iii), or it follows at once from (ii) ; for the least term 

— 2 f1 7 L\e~' a,n ( m 2m 

v w W f) 
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N m* 

-%/ 



agreeing with (iii). 

The following Table, to replace that in vol, xv. p. 432, was calculated in 
the way indicated above. 


X. 

Least term. 


Number of decimals 
directly obtainable 

1 

.. 3rd 


from the formula, 
2 

2 

.. 7 tli 


5 

5 

.. 16th 


13 

10 

32nd 


27 

20 

63rd 


54 

50 

157th 

.. 

136 

100 

314th 


273 

200 

629th 


546 

500 

1571st 


1365 

1000 

3142nd 


2729 


The number of decimal places practically obtainable is limited by the 
difficulty of calculating the Bernoulli’s numbers, of which only thirty-one 
have been hitherto obtained. By means of these, however, 156 decimals 
could be obtained of y when #= 100 ; and by deducing a few of the subse¬ 
quent terms, each from its predecessor (knowing their ratio), 20 places 
more could he obtained without difficulty. 

It is clear therefore that Mr. Shanks’s values of y obtained from #=500 
and #=1000 ought to agree beyond the 59th decimal, if correctly calcu¬ 
lated. The author at first supposed that the want of agreement was due 
to an insufficient number of the terms involving the Bernoulli’s numbers 
having been included, and he therefore undertook the calculation of this 
portion of the expression for #==500 and a?=1000 to 100 decimal places ; 
the results, however, still showed a difference in the 59th place. 

In order to determine where the error existed, the same portion of the 
constant was calculated also to 100 places, both from #=100 and #=200 ; 
all the calculations were performed wholly in duplicate, and so much care 
was taken that the author felt a strong conviction of their accuracy. It 
should be noticed that the agreement of all the four results would not 
necessarily prove the accuracy of that value of y; for any error made in 


the calculation of the harmonic series 1 -f- 


4-would merely pro- 
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duce the same error in 1+i . . . + ^-, l+~ * . . + 77 —* and 1+i . . . 

+ j q ' q q* supposing the calculation of the reciprocals beyond to 
have been accurately performed. 

It was therefore evident that the only means of ensuring freedom from 
error in the harmonic series was to recalculate it. The lowest value of x 
which would suffice for a verification was 100 ; the author therefore calcu¬ 
lated the value of the series l+~ • '• • 1° 100 places of decimals, 

2 I uu 


and the result was found to agree to that extent with that given by Mr. 

Shanks in the paper previously referred to; the value of 1 + ^ . . . +—* 

2 oU 


was also found to be correct. 

It may be mentioned that the calculation was abbreviated by a simple 
artifice, suggested by Oettinger (Crelle’s Journal, t. lx. p. 376), which will 
be easily understood from an example. Suppose the sums of the recipro¬ 
cals of the odd and even numbers up to 50 are known, and it is required 
to find the sum of the reciprocals up to 100. 

Let 


, ,1,1 
°- 1+ 3 + 5• 

1 

* ’ 49* 

/3=i+I + \ . 

1 

* ‘ 50 5 


then 


so that 


"+/ 3 =l+i + l +i+j . 

2 2 4 (5 98 100’ 


1+i+i . . . +—+— 

2 8 99 100 

= ^±£ + „+_L + _L. . . f±, 

2 51 53 99 


and the calculation of the reciprocals of the even numbers between 50 and 
100 is rendered unnecessary. 

The values of log 2 and log 5, which were required to form log 100, 
log 200, log 500, and log 1000 were taken from Mr. Shanks’s ‘Rectifi¬ 
cation of the Circle 5 (they are also given in the Proc. Roy. Soc. vol. vi. 
p. 397), and the summations of the harmonic series from Mr. Shanks’s 
paper in vol. xv. p. 431 of the Proc. Roy. Soc. The results were as 
follows:—- 

From x — 100 :— 


y = '57721 

56649 

01532 

86060 

65120 


90082 

40243 

10421 

59335 

93992 


35988 

05770 

42799 

90853 

75858 


22362 

13134 

84454 

84292 

91195 . . 

• (A). 
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From x = 200 :— 


y = *57721 

56649 

01532 

86060 

65120 


90082 

40243 

10421 

59335 

93995 


35988 

05770 

91390 

77706 

70283 


52607 

02838 

27101 

38332 

64367 . . 

■ • (B). 

From x = 500 :— 






y = -57721 

56649 

01532 

86060 

65120 ] 


90082 

40243 

10421 

59335 

93995 


35988 

05771 

53864 

75089 

05973 


20064 

75314 

10504 

07812 

29873 . . 

• (O). 

From x = 1000 :—* 






y = *57721 

56649 

01532 

86060 

65120 


90082 

40243 

10421 

59335 

93995 


35988 

05772 

02455 

61942 

00398 


50309 

65017 

53150 

61852 

03044 . . 

• (D). 


The terms involving B 39 , B 3l , B 22 , and B ]9 were the highest used in the 
calculation of (A), (B), (C), and (D) respectively. 

It will be seen that (A), (B), (C), and (D) differ in two respects : (i) the 
50th figure in (A) is 2, while in (B), (C), and (D) it is 5, and (ii) all four 
values are totally different after the 59th figure. 

The first discrepancy (i) pointed to an error in the summation of the 
harmonic series. As the author had verified Mr. Shanks’s value of 

1 +§•••+-[it was practically certain that (A) was the correct value. 

To place this beyond all doubt, however, the author calculated y from x 
= 50 to 57 places of decimals, and the result entirely confirmed (A). It 
follows, therefore, that Mr. Shanks has made an error of 3 in the 50th 

place in the calculation 200^ S0 we have the fol¬ 
lowing errata (vol. xv. p. 431) :—In 1 + TjT • • • +2^0’ ^ ie & rou P 

five should be . . . 53024 . . . instead of. . . 53027 . . . 

1 1 

In 1 +~2 ° • • + 5 QQ tenth group should be . . . 39677 . . . instead 
of. . . 39680 . . . 

In 1 . . . + the tenth group should be . . . 70880 . . . instead 

of. . . 70883 . . . 

With regard to the second error (ii), which causes the disagreement of 
all the figures after the 59th, it is clear that no error in the Bernoulli's 
numbers could produce such a discrepancy, as the terms involving the 
latter vary their position in each calculation, so that an error in any one 
of them could not affect the same decimal in (A), (B), (C), and (D). 
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By subtraction we find : 

(B) - (A) = . . . (60 ciphers) . 48590 86852 

94425 30244 89703 42646 

54039 73172 .'(E) 

(I)) _ (C) = . . . (60 ciphers) . . . 48590 86852 

94425 30244 89703 42646 

54039 73171 . .. ..(E) 


which only differ by a unit in the 100th decimal place. 

Leaving out of consideration the terms involving Bernoulli’s numbers, 
the difference between the series when # = 200 and when #=100 is 

(i + y • • • +200“log 200 ) — (l + f ‘ * ’ + 100 ~ log 10 °) 

= 1oI + T 02 ' ‘ ' + 200 “ los 2 ‘ 

Similarly, the difference between the series when # = 1000 and when 
x = 500 is 

1 i i 

~501 + 502 " ' + 1000 i0g 2- 

Now as it is extremely improbable that two errors, exactly equal in 
amount, should have been made in the calculations of 

1 1 i i i 1 

101 + 102 * ' * + 200 and 501+502 ' * * + 1000’ 

we have very strong evidence that the value of log 2 is inaccurate, and 
that (E) is the correction to be applied to it. 


By subtracting (C) from (A), we obtain 

(C) — (A) = . . . (59 ciphers) ... 1 11064 84235 

30114 97702 62179 26049 23519 

38678 .(F) 


The difference of the corresponding series (omitting as before the terms 
involving the Bernoulli’s numbers) 

i 1 1 _ 

■”101 + 102 # * * + 500 “ og 5 * 

Having only this one* difference-result involving log 5, it is impossible 
to decide from (A), (B), (C), and (D) whether the harmonic series or log 
5 or both are in error ; but the following reasoning places it beyond all doubt 
that (F) is a correction to log 5, and that the sum of the harmonic series is 
correct. 

* By subtracting (B) from (D) we might get another, but the portion 2^1+20^ 
+50O k armon i c series, as well as log 5, would he common to both.—dune 1C, 
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Mr. Shanks computed log 2 and log 5 from formulae of the form 

Log 2 = 2 (7 P + 5 Q + 3 R), 

Log 5 = 2 (16 P + 12 Q + 7 R), 

P, Q, and R being infinite series (Proc. Roy. Soc. vol. vi. p. 397; ‘ Rectifica¬ 
tion of the Circle/ p. 88). 

If, therefore, any error was made in the calculation of P say, it would 
produce errors in log 2 and log 5 proportional to 7 and 16. On trial it 
was found that sixteen times (E) was equal to seven times (F), the difference 
being only such as an error of a unit in the 100th decimal of (E) or (F) would 
produce. This afforded a moral proof that Mr. Shanks had made an 
error equal to one-fourteenth of (E) in the calculation of P, or 

1 1 1 
31 + 3.31 3 + 5.31 5 + • • • 

which has rendered his values of log 2 and log 5 incorrect, and that (with 
the exception of the error previously noticed) the harmonic series was 
summed correctly. 

Since log 3 was calculated from the formula 

Log 3 = 2 (11 P + 8 Q + 5 R), 

its value is also incorrect, as also is that of log 10 (log 2 +log 5) and the 
modulus (the reciprocal of log 10). 

The values of all these quantities are given to 205 decimal places in 
vol. vi. Proceedings of the Royal Society, p. 397; but all the figures after 
the 59tli decimal place are incorrect in each case. 

The correct values to 100 places are :— 


Log 2=*69814 

71805 

59945 

30941 

72321 

21458 

17656 

80755 

00134 

36025 

52541 

20680 

00949 

33936 

21969 

G9471 

56058 

63326 

99641 

86875 

Log 3 = 1-09861 

22886 

68109 

69139 

52452 

36922 

52570 

46474 

90557 

82274 

94517 

34694 

33363 

74942 

93218 

60896 

68736 

15754 

81373 

20888 

Log 5 = 1*60943 

79124 

34100 

37460 

07593 

33226 

18763 

95256 

01354 

26851 

77219 

12647 

89147 

41789 

87707. 

65776 

46301 

33878 

09317 

96108 

Log 10=2-30258 

50929 

94045 

68401 

79914 

54684 

36420 

76011 

01488 

62877 

29760 

33327 

90096 

75726 

09677 

35248 

02359 

97205 

08959 

82982 


It is to be observed that the above value of log 3 is not quite as well 
determined as the others ; the calculations in regard to Euler’s constant 
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form a real verification of log 2, log 5, and log 10 ; they also verify P, 
Q, and R ; but an error in log 3 in the transcription of P, Q, and R, or 
their multiplication by 11, 8, and 5, or in the final addition and multipli¬ 
cation by 2, would not be detected. 

In the above logarithms the last figure may be in error to the extent of 
one or two units. 

The value of Euler’s constant to 100 decimal places is— 


rH 

CM 

t- 

Ip 

II 

56649 

01532 

86060 

65120 

90082 

40243 

10421 

59335 

93992 

35988 

05767 

23488 

48677 

26777 

66467 

09369 

47063 

29174 

67497 


The last figure here also may be in error to the extent of one or two 
units. 

It will be observed that Mr. Shanks’s value* of y for x = 500 differs 
from (C) in the 65th decimal place, and that his value for x~ 1000 differs 
from (D) in the 73rd place. This is caused by an inaccurate value of B ]3 

8553103 

having been made use of. The correct value of B 13 is-^-; but Euler, 


who first calculated it, made it -g- 0 Acta Petropolitana ’ for 1781, 

p. 46), and this incorrect value is given in the *Penny Cyclopaedia’ 
(Article “ Numbers of Bernoulli ”) and probably in other places. 

The values of the first thirty-one Bernoulli’s numbers are given in a 
paper by Ohm (Crelle’s Journal, t. xx. p. 11), and B 13 is given correctly 
there. The agreement of the values of Euler’s constant contained in this 
paper (when the logarithms of 2 and 5 are corrected) afford a complete 
verification of the Bernoulli’s numbers as far as B 2e , and partial verifica¬ 
tions of the rest. 

The difficulty and inconvenience of making calculations to so many de¬ 
cimal places is sufficient to warrant the publication of the values of the 
positive and negative parts of the portion of the series involving the Ber¬ 
noulli’s numbers, in case any one should desire to repeat any part of the 
calculation. We have 


1 1 
r_1+ 2 + T 


, 1 IB, 

•+ir- lo s*-2i+2?- 


3.3 


and if m denote the sum of the terms of the same sign as the harmonic 
series, and n the sum of the terms of the same sign as the logarithms, 
viz. if 

B x b 3 b 5 

m 2 x 0 6x 6 10# 10 


1 b* 

= 2^+& + 8x + 


* Proc. Roy. Soc. vol. xv. p. 432. 
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then, when *= 100 : — 


w=*00000 

83333 

33337 

30158 

73773 

44885 

67821 

37321 

67823 

08773 

00082 

30639 

33761 

49846 

36254 

14224 
589 . . . 

82920 

15089 

52013 

00966 

n= ’00500 

00000 

83333 

33375 

00000 

21092 

80052 

53967 

77636 

59871 

21105 

34750 

25738 

21912 

47654 

32311 

133. . . 

34887 

61384 

97230 

71640 

When *=200 

*00000 

20833 

33333 

39533 

73016 

61283 

59538 

59711 

91078 

12258 

22890 

92137 

92329 

17052 

51095 

23490 

044 . . . 

89894 

78957 

35141 

06774 

n = *00250 

00000 

05208 

33333 

49609 

37505 

14960 

84887 

28877 

09510 

05685 

64006 

82170 

11770 

13023 

01074 
042 . . . 

35246 

38458 

75084 

75675 

When x — 500:— 

m = *00000 

03333 

33333 

33358 

73015 

87309 

34487 

73462 

42678 

21147 

87864 

83345 

53789 

61875 

16264 

59799 

196 . . . 

20761 

22352 

44796 

75809 

^ = •00100 

00000 

00133 

33333 

33344 

00000 

00008 

63960 

92825 

14227 

06303 

41363 

50645 

63522 

06215 

12108 

386 . . . 

84946 

96404 

61971 

02278 

When a-=1000 

m=-00000 

00833 

33333 

33333 

73015 

87301 

59487 

73448 

77428 

21067 

82137 

32165 

00876 

22580 

99464 

52882 

922 .. . 

37182 

41954 

63931 

48991 

n=*00050 

00000 

00008 

33333 

33333 

37500 

00000 

00210 

92796 

09323 

93526 

00040 

82093 

23718 

92820 

00915 

30332 

37713 

49061 

83360 


485 . . . 

The calculation was performed to 105 places, and the last two figures 
have been rejected. 
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Postscript. 

Received June 14, 1871. 

After the completion of the above paper, the author found that Mr. 
Shanks had, in second, third, and fourth supplementary papers on the 
Constant*, extended his calculations so as to determine y from #=2000, 
5000, and 10,000. 

The values so obtained all differ in the sixtieth decimal; in fact the 
higher x is taken, the further from the truth are the results, as the errors 
in the logarithms are multiplied by larger factors. 

The calculation for a?—2000 affords a verification of the error of log 2 ; 
for on subtracting the value of y (#=1000) from y (#=2000), we obtain 
(after correcting B 13 ) a result agreeing with E to 80 decimal places 
(which is as far as Mr. Shanks has calculated the latter value of y), with 
the exception of a difference of a unit in the 73rd figure—an error pro¬ 
bably in the summation of the harmonic series from 1000 to 2000. 

The values for #=5000 and #=10,000 are (besides the errors pre¬ 
viously noticed) inaccurate from the 62nd figure. 

VIII. “ Records of the Magnetic Observations at the Kew Observatory. 
No. IV.—Analysis of the principal Disturbances shown by the 
Horizontal and Vertical Force Magnetometers of the Kew Ob¬ 
servatory from 1859 to 1864.” By General Sir Edward Sa¬ 
bine, K.C.B., President. Received June 15, 1871. 

(Abstract.) 

This paper exhibits an analysis of the principal disturbances recorded 
by the horizontal and vertical force self-recording magnetometers of the 
Kew Observatory in the years 1859 to 1864, showing the progressive 
diminution in the number and value of the disturbances from a maximum 
in 1859 to a minimum in 1863, being the first moiety of the “ decennial 
period and showing also the distribution of the disturbances, increasing or 
diminishing the respective forces, in the several years, months, and hours. 

IX. <c Amended Rule for working out Sumner's Method of finding 
a Ship's Place.” By Prof. Sir William Thomson, F.R.S. 
Received June 15, 1871. 

In my previous communication on this subject (i anted , p. 259) I 
described a plan according to which, in the first place, two auxiliary 
lines were to be drawn on the chart, from two sets of numbers taken 
out of a proposed Table, and then Sumner’s line (the line on which the 
observation shows the ship to be) was to be interpolated, dividing the 
space between them in the proportion of the differences of the sun’s decli- 

* Proc. Hoy. Soc. vol. xvi. pp. 154, 299, vol. xviii. p. 49. 



